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1. Introduction 

Heavy quarks are related to some of the most exciting physics studies at hadron colliders. 
It is very likely that the third quark family has a special role in the breaking of electroweak 
symmetry. The mass of the top quark measured at the Tevatron [1] is a sensitive probe of 
the theory for the breaking of the electroweak symmetry, and it constrains, for example, 
together with other electroweak precision measurements the mass of the Higgs boson. 

So far, top quarks have only been produced at the Tevatron. Detailed studies of the 
properties of the top quark will be a main theme for the experiments at the Large Hadron 
Collider. The LHC is often termed a "top factory" since it is capable of producing many 
such particles per second. The top-pair production cross-sections will be measured with 
a negligible statistical uncertainty in comparison to the most optimistic predictions for 
the attainable accuracy of theoretical calculations. Systematic experimental uncertainties 
could be nevertheless sizeable. For example, the CMS collaboration anticipates to measure 
the top quark cross-section with an early systematic uncertainty of 10% to 20% depending 
on the decay channel of the top quarks [2]. These systematic errors may be further improved 
with a large integrated luminosity. 

A large body of work has been devoted to obtaining precise theoretical estimates for 
heavy quark cross-sections. Next-to-leading-order (NLO) QCD corrections for the spin- 
averaged cross-section have been computed in Refs. [3-5]. Calculations with the full spin- 
dependence of the heavy quark decays were performed in Refs. [6,7]. The effect of soft-gluon 
resummation at the leading and next-to-leading logarithmic approximations was accounted 
for in Refs. [8—11]. Soft gluon resummation effects beyond the next-to- leading logarithmic 
approximation were included in Refs. [12]. NLO QCD calculations and parton shower event 
generators have been matched using the MC@NLO approach in Ref. [13]. 

A recent theoretical estimate for the total cross-section at NLO in QCD at the LHC 
is [14]: 

4 LO (LHC,m top = 171GeV) = 875+^[;;;^] (scales)^gjg (PDFs) pb . 

Similar analyses have been performed in Refs. [15]. The above theoretical uncertainty is 
marginally sufficient for a comparison with the projected systematic experimental errors. 
It will be important to improve upon it by performing a complete NNLO calculation. 

The only two cross-sections for hadron collider processes which have been computed 
beyond NLO in QCD are Drell-Yan lepton-pair production [16-19] and Higgs boson pro- 
duction [17,20-26]. It is interesting that while in Drell-Yan production the NLO theoretical 
estimate from varying the renormalization scales turns out to be reliable, this is not the case 
in the gluon initiated process of Higgs boson production. It has also been observed that the 
NLO calculations for Higgs production fail for efficiencies when experimental cuts vetoing 
radiation are applied [27]. Similar cuts must be applied in various analyses (e.g. [28]) for 
top-production when this is a background process for other interesting signals. 

While there are many commons among top-pair production and Higgs boson produc- 
tion, such as the dominant contribution from gluonic initial states and a heavy invariant 
mass being produced in the final state, there are good signs that the top-pair cross-section 
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may have a better converging perturbative expansion; for example the NLO K-factor [29] 
is smaller for top production than Higgs production. It is also very encouraging that 
event generators with NLO matching such as MC@NLO are in very good agreement with 
NNLO [27] in estimating experimental efficiencies for Higgs boson production. Neverthe- 
less, it is necessary that any good intuition about higher order effects beyond NLO is 
verified explicitly with an NNLO calculation. 

There has been a significant recent effort towards computing the NNLO corrections. 
This is a computation which requires several ingredients. The virtual 2-loop QCD correc- 
tions in the limit s, \t\, \u\ S> rriq were computed in Refs. [30,31] using a factorization the- 
orem [32] and explicit two-loop computations for massless parton scattering [33]. The full 
two-loop amplitude for qq — > QQ was computed numerically in Ref. [34]. The correspond- 
ing analytic result for the sub-amplitude with fermionic one-loop insertions in propagators 
was computed in Ref. [35]. 

A lot more work is required for a complete NNLO cross-section. For example, the 
NNLO calculation should include the computation for the cross-section pp — > itjet + X 
at NLO, which has been recently computed in Ref. [36]. However, at NNLO there exist 
additional infrared singularities arising from one or two massless partons being unresolved. 
For the result of Ref. [36] to be utilized, a new subtraction formalism for NNLO calcula- 
tions is required, extending non-trivially the theoretical knowledge gained in the NNLO 
computations for e + e~ — > 3jets [37] and pp — > H + X [24]. An alternative approach would 
be to apply the methods followed in Refs. [19, 22, 38] used for Higgs boson and Drell-Yan 
production. In both approaches there are several new problems which need to be addressed, 
and success is by no means guaranteed. 

An important ingredient of the NNLO computation is the one-loop amplitudes for 
gg — » tig, gg — > ti and the likes replacing gluons with light quarks in the external states. 
In the absence of a subtraction method, these amplitudes need to be computed through 
C(e 2 ) in the dimensional regulator e = 2 — |. These amplitudes have been computed to 
sufficiently high order in e in Ref. [39] . The one- loop squared contribution for qq — > it was 
recently given in Ref. [40]. 

In this paper we compute the gluonic one-loop amplitude gg — > it using fully numer- 
ical and analytical methods. This is a rather humble computation in comparison to the 
remaining problems for obtaining a NNLO cross-section. However, many of the problems 
which are foreseen for later tasks, especially related to the numerical efficiency of comput- 
ing the one-loop contributions for sub-processes with five external legs in the presence of 
additional phase-space singularities, may benefit from an efficient evaluation of the loop 
amplitudes. 

Our goal in this paper is to achieve an efficient evaluation of the squared one-loop 
amplitude for gg — > QQ. Recently, powerful new methods [41-49] have emerged for the 
computation of one-loop amplitudes based on the groundbreaking work of Refs. [41,44]. 
These methods have started to make an impact in phenomenology [50] but are still lacking 
the maturity of other approaches [51-55] which lead to impressive applications of an ex- 
traordinary technical difficulty. Due to their conceptual simplicity, however, it is realistic 
to anticipate a significant impact with the new methods. 
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We have performed two independent computations of the one-loop gg — ► ti amplitude 
through order 0(e 2 ). First, mostly for testing purposes, we use the method of Davydychev 
[56, 57] which is analytic. Our main computational method is a new extension of the the 
method proposed in Ref. [58], which was also analytic. Here we show how to achieve a 
numerical implementation. 

Ref. [58] proposed to perform reductions of tensor integrals in two stages: first, a re- 
duction in exactly four dimensions where every propagator is modified by adding a common 
mass parameter, and second, a mapping of the results in the four dimensional reduction 
to a reduction in D = 4 — 2e dimensions. It was proposed then that the four dimensional 
reduction could be performed using the method of [41]. For the mapping to work, the 
analytic behavior in the mass parameter was required. Here, we find that a numerical re- 
duction for a fixed number of values of the mass parameter is sufficient to reconstruct the 
full dimensional dependence of the master integral coefficients at each phase-space point. 

A very elegant solution to the problem of a D-dimensional reduction for one-loop 
amplitudes directly and not just tensor integrals was given in Ref. [46]. In fact, the 
gg — > tig one-loop amplitude which is needed for the NNLO computation of top-pair 
cross-section was computed recently in Ref. [59]. The method of Ref. [46] capitalizes on 
unitarity ideas [60-65] and has the advantage of being able to use tree amplitudes for the 
reconstruction of master integral coefficients. 

Unitarity methods require the evaluation of one-loop amplitudes in renormalization 
schemes such as the four dimensional helicity scheme (FDH) [66] where external particles 
are in four dimensions. The FDH scheme has been successfully used for the evaluation 
of one-loop and two-loop amplitudes in QCD. However, it is known to be inconsistent for 
the two-loop amplitude gg — > h [67], and at higher loop orders the renormalization of the 
strong coupling is not equivalent to the conventional Dimensional Reduction scheme [68]. 
It is not yet clear to us whether the FDH scheme is a consistent scheme at NNLO for 
top-pair production. Our method allows us to compute the required one-loop amplitudes 
in conventional dimensional regularization (CDR) [69] by using a projection method in 
order to extract the spin dependence of the amplitude, similar to what it was used for in 



The paper is organized as follows: in section 2 we introduce our notation. In section 3 
we discuss the analytic method we used to check our results and our new numerical method. 
In section 4 we present our results for the gg — > QQ amplitude. Finally we conclude with 
section 5. 

2. Notation 

We consider the scattering process 



Ref. [70]. 



g{pr, ai) + g(p 2 ; a 2 ) + Q(j> 3 ; i 3 ) + Q(P4', U) -> 



(2.1) 



where we take all momenta pi incoming so that 



Pi + P2 + P3 + P4 = 0. 



(2.2) 
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Figure 1: Gluon amplitude for heavy quark production 



We denote by a\ , a-i = 1 . . . iV 2 — 1 and by i 3 , i$ = 1 . . . N the color of the external gluons 
and external quarks correspondingly. In addition, the momenta satisfy 

p 2 =p 2 = 0, P l=p\ = m 2 Q , (2.3) 

where we apply the approximation that the heavy quarks are on-shell and thq denotes the 
heavy quark mass. We define the Mandelstam variables, 

s=(pi+p 2 ) 2 , t=(p 2 +p?,) 2 and u = (p 1 +p 3 ) 2 ; (2.4) 

they satisfy 

s + t + u = 2m 2 Q . (2.5) 

We present here the unrenormalized QCD amplitude with Ni = light quark flavors 
and Nh = heavy flavors different than the flavor of the produced quark pair. In our main 
result, we also drop diagrams with self-energy corrections on the external legs, envisaging an 
on-shell renormalization procedure for the heavy quark mass and a decoupling of the heavy 
flavor from the renormalization group evolution of the strong coupling. The general result 
for the unrenromalized amplitude, having iV/,,JVj ^ will be given in Section 4. Infrared 
and ultraviolet singularities are regularized using conventional dimensional regularization 
performing both the 7-matrix algebra and loop integrations in D = 4 — 2e dimensions [69] . 

We write the amplitude as a perturbative series in the bare coupling a s , 

M gg ^ QQ = 4ira s 

The tree and one-loop amplitudes are expressed in terms of the color functions 

d=T%T%, c 2 = T%T%, C3 = 5 a ^5 Hl3 , (2.7) 

as 

L=l 

where n = 0, 1. The SU(N) generators in the adjoint representation are normalized as 
tr (T a T b ) = \5 ab . 



M 



(o) 



gg^QQ 



+ S r < 1 + E) 



m. 



M 



(i) 



99- 



+ 0{al). 



(2.6) 
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We present the amplitude for a specific choice of gluon polarizations, requiring that 
the polarization vector of one gluon is perpendicular to the momentum of the other gluon, 



e± -P2 = £2 -Pi = 0. (2.9) 

With this constraint on gluon polarizations the color coefficients Ml can be expressed in 
terms of a basis of "spin structures" as 

A4 B) =E$! r *' (2-10) 

i 

where 



ri 


= m Q ei- e 2 v(p 4 )u(p 3 ) 


r 2 


= m Q ei- p 3 e 2 • p 3 v(p 4 ) u(p 3 ) 


r 3 


= ei • e 2 v(p A ) p'iu(p 3 ) 


r 4 


= Cl • P3 £2 • £>3 «(P4) ^1 «(P3) 


r 5 


= £2 -P3V{P4) Au(p 3 ) 


r 6 


= £l -P3^(P4) ^fej) 


r 7 


= £2 -P3 v{p±) i\ A«(j?3) 


r 8 


= £l -J»3V(P4) <2 ?^1«(P3) 


r 9 


= mQv(p4) 4\ ^2u(p 3 ) 


Tio 


= mQv{pA) 4\ i>\ fou{jp 3 ). 



(2.11) 

In terms of these spin structures and the color basis of Eq. (2.7), the born amplitude 
is given by 

2 2 2 1 1 

iMf> = — r 3 ^r 3 + r r 6 T — r 10 

s u — m,Q u — rriq u — rriq tuq 

iMf ] = -r 3 + -^r 6 --^— r 10 , 

s t — rriq t — uig thq 
iM { 3 0) = 0. (2.12) 

In this paper we present the one-loop coefficients df 1 ^ through order 0(e 2 ), which we 
compute with both an analytic and a purely numerical method. The numerical method 
relies on extracting these coefficients by acting with a set of projectors on the one-loop 
amplitude. 

Let us multiply both sides of Eq. (2.10) by E- and sum over gluon and quark polar- 
izations 

E Mi n M = £ E $r*rt 

polarizations i polarizations 

7( n )n 



/ , a i,L L V> 



(2.13) 



where we defined the symmetric matrix to be 



polarizations 
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The coefficients di are then given by 



d) 



(n) 



(2.15) 



polarizations 



The matrix Tij and the inverse matrix T^ 1 are evaluated in CDR and we list their compo- 
nents in Appendix A. 

Once we know the coefficients di using Eq. (2.15), we can calculate the one-loop ampli- 
tude squared \M gg ^qQ\ 2 summed over external quark and gluon polarizations and averaged 
over initial state color and polarizations to be 



E 



|M 



(i) 



polarizations 



gg^QQ 1 



12 4 

^ — a Q 



256(1 - e) 2 



r(i + e) 



47T/T 



rn 



Q 



E 

L,M=1 



E E di > L ( e > ^ m< 3' tPx}) d j,M( e > m; "iq, {w}) r„ 



C LM (2.16) 



where we average over color and polarizations for the external gluons. We have also defined 
the color matrix Clm = tr (clcm), given by 



/ NC 2 F -\C F NC F \ 



Clm 



-\C F NCp NC F 
\NC F NC F AN J 



(2.17) 



with C F 



n z -i 

2N ■ 



We expand the one- loop amplitude squared in e, 

2<E 



polarizations 



(1) I 

gg^QQ 1 



4 T(l + e) 2 / 4ttm 2 



(l-ef 



rn 



Q 



A 4 \ + A 3 \ + A 2 \ + A 1 - + A 
e 4 e d e z e 



+0(e) 
(2.18) 

In Section 4 we present plots for A4, A3, A2, A± and Aq as a fuction of —t/s for fixed 
s = 16 m 2 -, and mg = 1. We also provide analytic results. 

When the heavy quarks are an intermediate state and a decay is considered, we simply 
need to replace the matrix with the analogous matrix for the specific decay in Eq. (2.16). 



3. The Method 

We have performed the calculation using two independent methods for the reduction to 
master integrals. Both methods are based on Feynman diagrams, which we generate using 
QGRAF [71]. Algebraic manipulations involving color and Lorentz indices are performed 
using FORM [72]. 
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3.1 Analytic approach 

In our first approach we perform an explicit analytic evaluation of all tensor integrals 
which appear in the one-loop amplitude. Following the method of Davydychev [56, 57] 
we write tensor integrals in terms of scalar integrals where the denominators are raised 
to integer powers of propagators and the dimensionality of the integral is increased by an 
even number. Schematically, 

k^ ... k^ 

- - \ 7 /'J-'-'/'' \ \ 

,;vi...v n 



I 



d°k- 



ill' 



T m v\...v n 

x f d D+2m k - 



m- 



l + Vi 



(3.1) 



where T^i—Vr are a \\ tensors of rank r that can be constructed from products of the 
metric tensor and the momenta m and Vj are non-negative integers and the coefficients 
;v\...v n are a ^ so integers. 
As a next step we perform a reduction of the dimensionally shifted integrals with extra 
powers of denominators to master integrals in the same dimension. These are typically 
scalar integrals with unit powers of propagators. 



2m 



k- 



IL (k + q 3 f- 



in' 



1 + 



I 



d D+2m k- 



IL 



(k + qj) 2 - m 2 



We perform this reduction using the program AIR [73]. 

Master integrals in D + 2m dimensions can be reduced in terms of master integrals in 
D dimensions [56,74,75]. 

l /■ : 



/ 



d D+2m k- 



IL 



(k + Qj) 2 - m 2 



[ d D k- 



m. 



l • 



(3.2) 



This concludes our analytic reduction of tensor integrals to master integrals. In Appendix B 
we present generic dimensional shift identities with arbitrary masses and arbitrary external 
momenta for box, triangle and bubble topologies, which we obtained with the method 
outlined in this section. 

This analytic method has worked fine for our application in this paper. However, we 
had to deal with rather large expressions in intermediate stages as well as in the final result. 
This can be neatly avoided with an alternative numerical method which we present here. 
We anticipate that this alternative approach will also facilitate the evaluation of the loop 
contributions to the NNLO cross-section from processes with five external legs. 

3.2 Numerical approach 

We consider the one-loop amplitude in D = 4 — 2e dimensions 

A(k) 



M(D 



f d D k 

2e) " J J^ft- 



(k + qi 



(3.3) 
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where k is the D dimensional loop-momentum. Following [62,63]. we may decompose the 
D dimensional integral into a 4 and aD-4 dimensional integral by writing k as 



k = I + n . 



(3.4) 



where I is the 4 dimensional part and fi is the D— 4 dimensional part for the loop momentum. 
We have, 

fc 2 = Z 2 -/A (3.5) 

In addition, we consider all combinations of external momenta qi in four dimensions, so 
that 

H-qi = 0. 

With this decomposition of the dimension and the loop-momentum, we can rewrite the 
integral in Eq. (3.3) as 



M(D = 4 - 2 C ) = *_ f (^p" 1 / d 4 l~ 



A(l,fi) 



(/ + ^) 2 -(M 2 + m 2 ) 



(3.6) 



In the above equation, the number of dimensions D = 4— 2e for the original loop momentum 
k is assumed to be bigger than four. Before we interpret the inner integration as an 
integration in four dimensions, the issue of the dimensionality of 7 matrices, as well as 
polarization vectors and spinors for the external states needs to be addressed. These are 
usually dealt with by performing the evaluation of the inner integrand using the FDH 
scheme [66]. As we explained in the introduction, we would like to avoid performing the 
computation in FDH. 

These issues do not arise when the amplitude is multiplied with the matrices Tj of 
Section 2 and the trace is taken. This is sufficient in order to reconstruct the full amplitude 
once the matrix Tij = tr (v^Tj^ and its inverse are known (appendix A). We then write, 



tr 



T\M{D = A-2e) 



vr 2 r(-e) 



d» 2 (^ 2 ) 



-e-l 



d 4 l- 



tr 




CDR 




( 



(3.7) 

where we have performed the trace in CDR before decomposing the loop momentum. 
Therefore, the traced integrand depends explicitly on e from the 7-matrix algebra in CDR. 



tr 



r\M{D = 4-2e) 



7r 2 r(-e) 



di? 



-e-1 



dH- 



CDR 



(3.8) 



In what follows we perform a separate reduction for each one of the CDR e coefficients 
in the numerator. The dependence on e due to the loop integration can only be made 
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manifest if the integration over /j 2 is performed. However, following the spirit of [58], we 
will avoid this integration. 

We now focus in the four dimensional inner integration. For every one-loop amplitude 
there exists a reduction 



dH- 



tr 




CDR 









^f4Gu 2 )Master£= 4 (fj 2 ) , 



(3.9) 



in terms of master integrals Master^ in the same dimension D = 4. The master integrals 
in four dimensions are the scalar tadpole, bubble, triangle, and box integrals with unit 
numerator in their integrand. One can enhance however, the basis of master integrals by 
including the scalar pentagon. With this enhanced basis the coefficients of the master 
integrals are simple terminating polynomials in /j, 2 [58,76]: 



^Gu 2 ) 



4 0) 



+ /i 2 fl[ 1) + ... + (rt H 



(3.10) 



The power max can be determined on simple dimensional grounds for each process and 
theory. We shall comment later on the methods that may be used to achieve this reduction. 
What is important to note now is that it can be performed numerically, since no divergence 
emerges. We then choose max different values of fi 2 at each phase-space point and form 
a system of equations from Eq. (3.10) 1 . We solve the system numerically, and obtain the 



values of the coefficients fij^ . 
>(0 



With the Q k already evaluated it is easy to obtain the reduction in D 
dimensions. We write, 



2c 



/ 



d D l- 



tr 




CDR 


e" 






H 2 + m 2 ^ 



7r 2 r(-e) 



dlJ, 2 (fj 2 ) 



dH- 



tr 




CDR 








fj, 2 + mfj 



J>« x -^L_ J V ( ^ 2) — Master^ 4 (, 2 ) 



= ^ x Master^ 



D=4-2e+2i / 2 _ 



(, 2 = o) • 



(3.11) 



In the last line we have achieved a reduction in terms of master integrals in 4 — 2e + 2i 
dimensions, where neither the master integrals nor the coefficients depend on the mass 
parameter fi 2 . These dimensionally shifted master integrals can be computed easily in 
terms of master integrals in D = 4 — 2e using the procedure described in Eq. (3.2). 

We now comment on the methods which we have employed to perform the four di- 
mensional reduction of Eq. (3.9). Our application is simple enough so that an analytical 



X A similar inversion technique was previously used for a different purpose in Ref. [46] 
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or numerical reduction using the program AIR [73] was possible. It is worth noting that 
the four dimensional system of integration by parts identities which we need to solve in 
Eq. (3.10) is much simpler than what it is required for an equivalent direct reduction in 
D-dimensions. The two reductions involve the same number of symbolic parameters but 
they differ in complexity ( we have introduced a new mass parameter /j 2 in Eq. (3.9) but at 
the same time we take e = 0). This is mainly due to the fact that the reduction coefficients 
are a simple polynomial in /i 2 in the four dimensional reduction, while in the D-dimensional 
reduction the coefficients are in general rational functions. 

The best reduction method which we have used in Eq. (3.9) is the one introduced by 
Ossola, Papadopoulos and Pittau [41]. This method finds the coefficients of master integrals 
by evaluating the loop integrand at a finite number values of the loop momentum. This is 
sufficient, since the integrand in four dimensions has a known functional form as a sum of 
a small number of rational functions in the loop momentum [41] . 

A different method than ours [46] for a D-dimensional reduction also exploits the known 
functional form of one-loop integrands [41] and it computes the coefficients of the master 
integrals purely from tree amplitudes [45]. Recently, the one- loop amplitude for gg — > tig 
was computed with this alternative method [59]. However, this method was used only in 
the FDH and the 'tHooft-Veltman schemes which we have decided to cautiously avoid as 
a first step. Finally, we note that our method is very close to what was described in [58], 
although an analytic approach was envisaged there. In [58], it was suggested that spinor 
integration [77] could also be used instead of the method of [41] for the four dimensional 
part of the reduction. The spinor integration method has developed significantly in the 
last couple of years [76,78-81]. However, it remains an outstanding problem to compute 
the coefficient of the tadpole master integral which vanishes when two propagators are cut. 

4. Results 

As we discussed in previous sections, we calculated the NLO gg — > QQ amplitude both 
analytically and fully numerically. In this section we are going to present both results for 
the unrenormalized amplitude. 

For our analytic results, we provide a Mathematica file, mastercoeff.m. We present the 
coefficients for the master integrals defined below for the unrenormalized amplitude with 
Ni = 0, Nh = and N = 3. The master integrals that appear in the amplitude are defined 
as 
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T _R UK 2 ^_ f d k 1 1 

r_ r r-o/' 2 \ _ f dPk 1 1 1 

7 " Tna2(m «' t} - y ^ f (fc+p 3 )2- m ^ (fc+^+p^.^ 

2 f ft 1 1 1 

8 " Tna2(m «' U) ~ J ¥ ( k +p 3 ) 2 -m 2 Q (k+ Pl +p 3 ) 2 -m 2 Q 



ig = Tria3(mg, s 



dPk 



m D l 2 k 2 — uIq (k + p\) 2 — nig (k + p\ + P2) 2 — "»g 
2 n_ [ d D k I 1 1 1 

i 10 -Boxl(m Q , S ,t)-y -D72 ^2 (A; + pi) 2 (fc+;pi+p2) 2 (A . +J , 1+p2+j , 3 )2_ m 2 



In = Boxl (ttt-q, s, u) 



ft 11 1 1 



m D / 2 k 2 (k + pi) 2 {k + p\+ p 2 ) 2 (k - p 3 ) 2 - m 2 Q 

I - Box2(m 2 s t) - [ 1 1 1 1 

12 Q ' ' y i7r D / 2 fc 2 (fe + pi) 2 (A; + pi + P3) 2 - fn 2 Q (k + pi + p 2 + P3) 2 ~ m 2 Q 

,0 \ f dp* /c 1 1 1 1 

il 3 = Box2(m , S,u) = / FTaTT^T^ v>71 ^ tt, ^ T 

v j m D/z fcz ^ + P2 y (fc + p 2 + p 3 ) 2 - (fe +pi +p 2 +P3) - m Q 

/d D k 1 
i^fc 2 ^ 



1 1 1 



'15 



Box3 



m 2 Q (k + pi) 2 - m 2 Q (k + pi+ p 2 ) 2 - m 2 Q (k + pi + p 2 + P3) 2 

2 _ /" ft 1 1 1 1 

(m Q ,s,u)- J k2 _ m 2 Q{k + pi y _ m 2 Q ( k+pi+ V2 y - m 2 Q (k- p 3 ) 2 

(4.1) 



We used the results in Ref. [82] for numerical values of the master integrals defined above. 
In order to define the coefficients B(i, L,j, k) in the file we provide, we write 

15 4 

d S = E E B(i,L,j,k)I j( *, (4.2) 
j=i fc=_i 

where dj 1 ^ are defined in Eq. (2.10), Ij are the master integrals defined above and k is the 
power of the dimensional regulator e. We give only the B(i,L,j,k) needed for the C(e 2 ) 
expansion of the amplitude. 

Our numerical results are presented in Fig. 2, where the coefficients Ai in Eq. (2.18) 
are plotted as a function of —t/s. In these plots, we have normalized uiq = 1 and we 
chose s = 16 itiq. These numerical results agree of course with our independent analytic 
evaluation. We have checked that the infrared poles of our amplitude agree with the 
universal infrared pole structure as predicted in Ref. [83]. The finite part of the one-loop 
amplitude agrees with the result of [39]. 

The results we present in Fig. 2 and in the Mathematica file correspond to zero num- 
ber Ni of light quarks and zero number of heavy quarks with a different mass than 
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the external quark. Contributions from the diagrams shown in Fig. 3 need to be added 
separately. These extra terms read: 




r(l),JVl _ 
gg^QQ 



r(l) 

gg^QQ,i 




+ 0(a 3 s ) 



(4.3) 



where 



iM 



iM 



(i) 



gg^QQ 



QQ 




12 Tadp (m 2 ) (1 - e) 2 + 2 mf s 



Bubl(s)r 3 



Tria3 (m 2 ,s)(2e - 3) 



-12m 2 Bub3 (m 2 ,s)(l - e) - sBub3 (mf,s)e T 3 . 



(4.4) 



where the definitions for the master integrals are given in Eq. (4.1). Note that in the above 
equations, rrii are the masses of the quarks for the internal loops and tuq is the mass of 
the external heavy quarks. 

5. Conclusions 

In this paper we compute the one- loop amplitude for gg — > QQ through order (e 2 ) and 
its square through order (e°). This is a contribution, albeit an easy one to compute, 
to the NNLO cross-section. We performed the computation in two different ways. We 
have developed here a new reduction method for one-loop amplitudes which yields the full 
dependence of master integral coefficients in the dimensional parameter in conventional 
dimensional regularization. We are currently applying this method for computing other 
one-loop contributions at NNLO from processes with five external legs. 
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Figure 3: (a)-(c) : Extra diagrams with internal massless quark loops, (d)-(f) : Extra diagrams 
with internal quark loops for quarks with mass rm. 
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A. Tij and T" 1 

In this appendix, we are going to give the elements of the matrices Tij, defined in Eq. (2.14) 
and its inverse r^- . Note that both matrices are symmetric. 



T n = 2(-l + e)(s + u + t)(s + 2u + 2t) 

T 12 = -—(a + u + t)(s + 2u + 2t) 
4s 

x (s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 
r i3 = 2(-l + e)(u - t){s + u + t) 

r i4 = — -(u - t)(s + u + t)(s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 
As 

Ti 5 = --(s + u + t)(s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 
s 

Ti 6 = --(s + u + t)(s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 
s 

r i7 = -(s + u + t)(s 2 + 2s-u + u 2 + 2st - 2ut + t 2 ) 
ris = tO + « + *)(s 2 + 2s-u + -u 2 + 2st - 2nt + t 2 ) 
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rig = 2(-l + e)(s + u + t)(s + 2u + 2t) 
ri 10 = -2(-l + e)(u - t)(s + u + t) 

r 22 = -^( s + u + + 2« + 2t) 

x (s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 2 

r 23 = - — (u-t)(s + u + t) 
4s 

x (s 2 + 2su + u 2 + 2st - 2ui + i 2 ) 

T 24 = -j^(u-t)(s + u + t) 

x (s 2 + 2su + u 2 + 2st - 2ui + t 2 ) 2 

r 25 = — + u + t)(s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 2 
4s^ 

r 26 = --r^(s + u + t)(s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 2 
4s^ 

r 27 = (s + u + t)(s 2 + 2su + u 2 + 2st - 2nt + t 2 ) 2 

16s 

r 28 = — (a + u + t)(s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 2 

T 29 = --^(s + u + t)(s + 2u + 2t) 

x (s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 

r 2 i = — (u-t)(s + u + t)(s 2 + 2su + u 2 + 2st-2ut + t 2 ) 
4s 

T 3 3 = -(-1 + e)(s + u - t)(s - u + t) 
T 34 = —(s + u-t)(s-u + t) 

x (s 2 + 2su + u 2 + 2st - 2ui + i 2 ) 
r 35 = - — ((«- t)(s 2 + 2su + u 2 + 2st - 2ut + t 2 )) 

T 36 = — -({u- t)(s 2 + 2su + u 2 + 2st - 2ui + t 2 )) 

r 37 = o 
r 38 = o 

T 39 = 2(-l + e)(u-t)(s + u + t) 
r 3 io = (-l + e)(s + u-t)(s-u + t) 

T 4 4 = (s + « - - « + t) 

x (s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 2 

r 45 = -g^a ((« - *)(« 2 + 2su + + 2st - 2ui + t 2 ) 2 ) 

T 46 = — — ((« - t)(s 2 + 2su + u 2 + 2st - 2ui + t 2 ) 2 ) 
8s 2 
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r 47 = o 
r 48 = o 

r 49 = - — ((u - t)(s + u + t)(s 2 + 2su + u 2 + 2st - 2ui + t 2 )) 

T410 = - — (s + u - t)(s - u + t) 

x (s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 
r 55 = -L(-3s 2 + 2es 2 - 2su - u 2 - 2st + 2ut - t 2 ) 
x (s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 

r 5 e = -J-2 (* 2 + 2su + « 2 + 2st - 2irf + * 2 ) 
x (2s 2 + 2su + + 2st - 2nt + t 2 ) 

r 57 = -^(-1 + e)(s + u + t)(s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 

r 58 = -(s + it + t)(s 2 + 2sit + it 2 + 2st - 2nt + t 2 ) 

r 59 = --((« + n + t)(s 2 + 2siz + it 2 + 2st - 2ui + i 2 )) 

r 5 i = -!-((s-2es + u-t)(s 2 + 2su + ti 2 + 2st-2nt + t 2 )) 

T 66 = ^2 (-3s 2 + 2es 2 - 2su - u 2 - 2st + 2ut - t 2 ) 

x (s 2 + 2su + u 2 + 2st - 2irf + t 2 ) 
r 67 = j(s + it + t)(s 2 + 2sit + u 2 + 2st - 2ui + t 2 ) 

r 68 = --(-1 + e)(s + u + i)(s 2 + 2su + u 2 + 2st - 2ut + i 2 ) 

r 69 = ~{{s + u + t)(s 2 + 2sit + it 2 + 2st - 2ut + t 2 )) 

r 6 i = -^-((-s + 2es + u-t)(s 2 + 2su + ti 2 + 2st-2ut + t 2 )) 

T 77 = — (-1 + e)(s + u - t)(s - u + t) 



x (s + u + t)(s 2 + 2su + u 2 + 2st - 2irf + t 2 ) 
-u + t) 

(s + u + t)(s 2 + 2su + u 2 + 2st - 2irf + t 2 ) 

T 79 = — -(s + u + - it + 2eu + t - 2et) 

4s 

x (s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 



F 78 = -7tH s + u ~ ~ u + *) 
lbs 



r 7 i = o 

r 88 = g^(-l + e){s + u-t)(s-u + t) 



- 16 - 



x (s + u + t)(s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 
-^(s + u + - u + 2eu + t - 2et) 
x (s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 



--2(s + u + i) 
s 

x (s 2 - 2es 2 + 2eV + 2su - Aesu + Ae 2 su - u 2 

+ 2eu 2 + 2st - Aest + Ae 2 st + 2ut - Atut - t 2 + 2et 2 ) 
4(-l + e) 2 (u-t){s + u + t) 
2(-l + e) 2 (s + u - t)(s - it + t) 



(s + u — t)(s — u + t) 
2es(s + u + t)(s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 

4(s + u-i)(s-u + i) 
~ (-1 + 2e)(s + it + t)(s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 2 
u — t 

es(s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 
8{u-t) 

~ (-1 + 2e)(s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 2 





2(u-t) 

e(-l + 2e)s(s + u + t)(s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 

2{u-t) 

~ e(-l + 2e)s(s + u + t){s 2 + 2su + u 2 + 2st - 2ttt + t 2 ) 

(s + u — t)(s — u + t) 
2e(-l + 2e)s(s + u + i)(s 2 + 2su + u 2 + 2st - 2ui + i 2 ) 
u - t 

~ e(-l + 2e)s(s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 
32(-2 + e)s(s + u-t)(s -u + t) 
~ (-1 + 2e)(s + u + t)(s 2 + 2su + u 2 + 2st - 2nt + t 2 ) 3 
8(u-t) 

~ (-1 + 2e)(s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 2 

64(-2 + e)s(u - t) 
~ (-1 + 2e)(s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 3 



16s 

(-1 + 2e)(s + u + t){s 2 + 2su + u 2 + 2st - 2itt + t 2 ) 2 



r 



28 



16s 

(-1 + 2e)(s + u + t)(s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 2 



r£ = 

r^io = o 

i s + 2u + 2t 

33 " es(s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 

i 8(s + 2n + 2t) 

34 ~ (-1 + 2e)(s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 2 



r 



i 



-i 

" e(-l + 2e)(s 2 + 2su + u 2 + 2st - 2nt + t 2 ) 

1 



r- 1 = 

36 e(-l + 2e)(s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 



r 



l 

37 



4 



r 



-i 



e(-l + 2e)s(s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 
4 



.38 



e(-l + 2e)s(s 2 + 2su + u 2 + 2st - 2ui + t 2 ) 
ti — £ 



r- 1 



e(-l + 2e)s(s 2 + 2su + u 2 + 2st - 2ut + i 2 ) 

r _i s + 2u + 2t 

3 10 ~ e(-l + 2e)s(s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 

! _ 64(-s 2 + es 2 - 2su + 2esu + u 2 - 2st + 2est - 2ut + t 2 ) 

44 ~ (-1 + 2e)(s 2 + 2su + u 2 + 2st - 2ut + i 2 ) 3 
p-i = g(tf - t) 

45 " (-l + 2e)(s 2 + 2su + n 2 + 2st-2ut + i 2 ) 2 
r -i = gfo - t) 

46 " (-l + 2e)(s 2 + 2su + n 2 + 2st-2ut + t 2 ) 2 

rjr 1 = o 
r 4 "8 = o 
r 4 "9 = o 
rii, = o 

r _i _ (s + u — t)(s — u + t) 

55 ~ ~ e(s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 2 
p_i _ (s + u — t)(s-u + t) 

56 " ~ e(-l + 2e)(s 2 + 2sit + u 2 + 2st - 2ui + t 2 ) 2 



r 



57 



4s 



e(s 2 + 2su + u 2 + 2st - 2ui + i 2 ) 2 



P-i = l£ 

58 e(-l + 2e)(s 2 + 2su + u 2 + 2st - 2nt + t 2 ) 2 

r 5 "9 = o 



r 



510 



1 



e(-l + 2e)(s 2 + 2su + u 2 + 2st - 2ui + t 2 ) 
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r _i _ (s + u — t)(s — u + i) 



66 



1 67 



r -l 
1 68 



1 610 



77 

-1 



^78 



p-1 
1 79 



p-1 
1 710 



p-1 



-^89 



p-1 
1 810 



-^99 



p-1 
1 910 



e(s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 2 
As 

' e(-l + 2e)(s 2 + 2su + u 2 + 2st - 2nt + i 2 ) 2 
4s 

' e(s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 2 
1 



e(-l + 2e)(s 2 + 2su + u 2 + 2st - 2nt + £ 2 ) 
! _ 8(-2s 2 + 2es 2 - 2su - u 2 - 2st + 2ut - t 2 ) 



<- 


-1 + 2e)s(s + u + t)(s 2 + 2su + u 2 + 2st - 


2ut + t 2 ) 2 




8(2s 2 + 2su + u 2 + 2st - 2ui + t 2 ) 




<- 


-1 + 2e)s(s + u + t)(s 2 + 2su + u 2 + 2st - 


2ut + t 2 ) 2 




2(u - t) 




<- 


-1 + 2e)s(s + u + t)(s 2 + 2su + u 2 + 2st - 
4 


2ut + 1 2 ) 


<- 


-1 + 2e)s(s 2 + 2su + u 2 + 2st - 2tti + t 2 ) 






8(-2s 2 + 2es 2 - 2su - u 2 - 2st + 2ut - 


~t 2 ) 



5 e(-l + 2e)s(s + it + t)(s 2 + 2su + u 2 + 2st - 2ut + i 2 ) 2 
! 2{u-t) 



e(-l + 2e)s(s + u + t)(s 2 + 2su + u 2 + 2st - 2trf + i 2 ) 
4 



e(-l + 2e)s(s 2 + 2su + u 2 + 2st - 2ui + i 2 ) 

! (S + U — t)(s — U + t) 



2e(-l + 2e)s(s + u + t)(s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 
u — t 

e(-l + 2e)s(s 2 + 2su + u 2 + 2st - 2ut + t 2 ) 



r -i = s + 2u + 2t 

10 10 e(-l + 2e)s(s 2 + 2su + u 2 + 2st -2ut + t 2 )' { ' ' 



B. Generic Dimensional Shift Identities 

In this appendix we are going to give generic dimensional shift identities for box, triangle 
and bubble topologies with arbitrary masses for different propagators. For this purpose 
let's define the following scalar integrals 



/d D k 1 1 

-W [{k + m ) 2 - m 2 ] [{k + qi ) 2 - m f 



1 1 



P + (?2 )2_ m 2] [{k + q3) 2_ m 2y 
[ d ° k 1 1 1 

Tri{m ,m 1 ,m 2 ,qo,qi,q2;-L>) = / . n/ 2 777— To 2TT71 - vi 2TT77 - Y2 21 > 

7 ?7r L V i [(fe + go) 2 - mg] [(fe + gi) 2 - mf] [(fe + g 2 ) 2 - m|] 
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Bub(m ,mi,g ,gi;-D) = J 
Tadp(m;£>) = J 



1 



1 



[(fe + g ) 2 - ml] [{k + gi) 2 - m\] 
d D k 1 



ivr D / 2 fc 2 



m, 



2 ' 



where in the definitions of the scalar integrals, the order of the arguments correspond one 
to one to the order of the propagators. 
For the bubble topology we have 

Bub(m o ,mi,0, qi;D + 2) 



x ^Bub(mo, mi, 0, gi; D) 



2(D-l)qj 
ml + (ml - g 2 ) 2 - 2m 2 (m 2 + g 2 ) 



+m\ 



Tadp(mi; D) + Tadp(m ; D) 



Tadp(mi; D) + Tadp(mo; D) 



m n 



Tadp(mi; D) + Tadp(mo; D) 



(B.2) 



For the triangle topology 
Tri(m ,mi,m 2 ,0,gi,g 2 ;-D + 2) 



x|Bub(mi,m 2 ,gi,g 2 ;-D) 



( J D-2)[g4 + (g 2 -g 2 2 )2_2g 2 (g 2 + g 2 2 ) 
(m 2 - ml)(q\ - q%) + {-2ml + m l + m l + Qi + <? 2 )<7i 2 ~ I12 



+Bub(m , m 2 , 0, g 2 ; D) 



+Bub(m , mi, 0, qv,D) 



m 2(<li +Q2~ Q12) + <l2{- 2m i + I1-Q2 + Q12) + m o(-1i +Q2 + Q12) 



m\(q\ + qj- q\ 2 ) + ™o(<7i -Q2+ Q12) + 9i(-2m 2 - q{ + g 2 + gi 2 ) 



2 „2 , „2 , 2 



4 2, 42 

m!g 2 + m g 12 



-2 Tri(m , mi , m 2 , 0, gi , g 2 ; D) 

(m\ - m 2 2 ){ql - gf) - (m\ + m 2 2 + q\ + q\)q\ 2 + q\ 2 



+ml 



m t + Q2Q12 + m\\q{ - g 2 - gi 2 ) - m 2 m 2 (gi + g 2 - gi 2 ) + q 2 {q{ - g 2 + gi 2 ) 



(B.l) 



where we defined 



Qij = ilj -ft) 2 



(B.3) 



(B.4) 
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For the box topology, the most general expression is quite lengthy. Therefore we will give 
the answer for q\ = 



Bax(mo,mi,m 2 ,m 3 ,0,qi,q 2 ,q 3 ;D + 2) = 



2(D - 3) \q\ 2 q 2 + ql 2 {q 2 l3 {-q 2 2 + qfa) - (qf 3 + q\ + 9|,)?3 + + g£(gf 3 + g| 3 «3 + TvM ~ T 2 3 ~ %)) 
x ^Tri(mo,m 1 ,m 3 ,0,q 1 ,q 3 ; D) - m\{q\ 2 - q 2 2 ){q\ 3 - q\) + m\{q\ 3 - q\f - (qj 3 - q 2 ){q\ 3 ql - q\ 2 ql) 

+ ml{ql 3 {q 2 l2 - q\ 3 - 2q\ + q\ 3 ) + {q\ 2 + q\ 3 - qj 3 )q 3 ) + m\ (ql(-2qf 2 + q\ + <?f 3 - gf ) 

+ 913(92-^3 + ^1)) 



+ TY\{mi,m 2 ,m 3 ,qi,q 2 ,q 3 ;D) 



2 2 2 42 222,222 
9l29l 3 9 2 - 9l 3 9 2 - 2gi 2 ^i 3 g 23 + 9l 3 9 2 923 



+ ™o(9r2 + (?13 - ?23) 2 - 2 9l2(^13 + Qh)) ~ Q12Q3 + ?129l3?3 + 9?2<723<?| 
+ m 2(9l3(9?2 - Ql3 ~ 2 Q2 + (Ql2 + Ql3 ~ Q23)4) 

+ »™§(-gf 2 + ^1(^13 - 923) + <7l2(<7?3 + 92 + Q23 ~ 2< &)) 



+ "h(9l 3 (9 2 + 923 - 9 3 ) + 923(92 - 923 + 9 3 ) + 9l 2 (-92 + 923 + 9 3 )) 

+ Tri(m ,m 1 ,m 2 ,0,,i, g2 ; J D) m 2 (q 2 2 ~ llf ~ m 2 2 (q 2 12 - q 2 2 )(q 2 13 - q 2 ) 
~ (Q12 ~ 9l)(-9? 3 9i + 9? 2 9i) + ™o(-9r2 + 9^(9?3 ~ Q23) 

+ 9? 2 (9? 3 + 92 + 9I3 - 2( &)) + m i(Qi2(Q2 ~ Q23 + <&) + (&{- 2( li3 -Q2 + Q23 + 9l)) 



+ Tri(m , m 2 , m 3 , 0, g 2 , 93; £>) 

+ Q2(Ql3- 2( l23) 



9| 



(91 - 9I5) 9?3( m 2 - 92) + m l(<l2 ~ Q23) 

1 2 m\{ql + q 23 ) + qf 2 (q 2 + q 23 ) + m 2 2 {-2ql 2 + g^ 3 + q\ + g| 3 ) 
- (-m? + + ^ 2 )g| + q\ 3 {q\ + g£ 3 - gf) + q\ 3 {q\ - q\ 3 + q\) + ^(-tff + 9^3 + 4) 

+ m i(9?2(9i - 923 + 9!) + d(- 2 9i3 -92 + 9I3 + 9i)) 



+ Box(m , m 1 ,m 2 , m 3 ,0, q±,q 2 , q 3 ; D) 



m 3(.Qi2 ~ qlY + m o 9l2 + (9r3 - 923^ - 2 9^(^3 + 923) 



-2n2 



2 \2 



+ "i2(9?3 - & + (9^392 - 9r293) 2 + fn\ q 2 + {q 23 - q 3 Y - 2^(^ 3 + gg) 



,2\2 



^2 2 2 2\2 



,2\2 



,2/2 



+ 2m? 



9l392(-92 + 923 + 93) + 93(- 2 92923 + 912(92 + 923 - 93)) 
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+ - (qf 3 - q1)(qf 3 qi - q 2 12 <&) + "if (g§(-2qf? 2 + q\ + g£ 3 - g§) + <7? 3 (<?1 - g| 3 + q\)) 

+ 2mg qj 3 q%(-qj 3 + <?1 3 ) - g^i + ^2(^3(^2 ~ 9?3 ~ 2 ^2 + ?2s) + (?12 + 9l3 ~ Sas)?!) 

+ m i(-^12 + ^1(^13 - 923) + 9l2(9?3 + 92 + ^23 ~ 2 li)) + l^i^l + Vizi^l ~ 2 ?23 + 83)) 
2 (9l3(?2 + 923 - 93) + ^23(^2 - 923 + 9s) + 9? 2 (-?2 + ?23 H 
" ™l(9?2 - 9i)(9?3 - 9l) - (9l2 - 32)(-?13«2 + 9l2^i) 



+ m 
+ 2m 3 



+ m 



1(912(92 - 9I3 + 9l) + g£(-2gi 3 - 9I + 9I3 + 9s)) 



(B.5) 
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